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SUMMARY 


Turbulent -heat -transfer  measurements  were  obtained  through  use  of  a 
•three-dimensional  axially  symmetric  nozzle  which  consists  of  an  inner 
shaped  plug  and  an  outer  cylindrical  sleeve.  Measurements  were  taken 
along  the  outer  sleeve  and  gave  flat-plate  results  that  are  free  from 
wall  interference  and  corner  effects  for  a  Mach  number  of  3*03  and  for 

a  Reynolds  number  range  of  5-6  x  10^  to  6.5  X  10^.  The  heat -transfer- 
coefficient  results  are  in  good  agreement  with  theoretical  analyses  and 
the  recovery-factor  results  are  in  good  agreement  with  extrapolations  of 
lower  Reynolds  number  data. 


INTRODUCTION 


The  design  of  supersonic  aircraft  and  missiles  requires  engineering 
information  about  heat-transfer  coefficients  and  recovery  factors  that 
extend  over  a  wide  range  of  Reynolds  numbers.  Each  of  the  existing  experi¬ 
mental  techniques  for  studying  heat -transfer  characteristics  appears  limited 
in  various  ways.  A  test  setup  which  utilizes  a  flat  plate  for  the  testing 
surface  is  subject  to  corner  and  edge  effects  and  also  shock  interference 
if  high  Reynolds  numbers  are  obtained  by  extending  the  plate. 

One  existing  technique  which  avoids  these  effects  is  that  of  the 
hollow  cylinder  mounted  centrally  in  the  test  section.  However,  the 
apparatus  is  subject  to  the  tube  choking  at  low  Mach  numbers  and  deviation 
from  flat-plate  conditions  (the  boundary- layer  thickness  approaching  a 
sizable  fraction  of  the  cylinder  radius)  when  the  Reynolds  numbers  are 
low.  In  addition,  there  are  differences  in  the  heat -transfer  character¬ 
istics  for  the  inside  and  outside  surfaces  of  the  cylinders. 

The  purpose  of  this  paper  was  to  present  results  obtained  by  using 
a  technique  that  would  avoid  the  above-mentioned  effects  and  permit  sim¬ 
ple  and  accurate  measurements  of  local  heat-transfer  coefficients  and 
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recovery  factors  over  an  extended  range  of  both  Reynolds  and  Mach  numbers. 
The  development  of  this  technique  for  Mach  number  3-03  and  the  results  of 
preliminary  heat -transfer  measurements  for  the  cooling  of  the  surface  by 
the  airstream  are  the  subject  of  this  paper. 
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SYMBOLS 


specific  heat  of  the  sleeve  material,  Btu/lb-°R 
specific  heat  of  the  air  at  constant  pressure,  Btu/lb-°R 


acceleration  due  to  gravity,  ft/sec^ 
heat-transfer  coefficient,  Btu/ft^-sec-°R 
heat  conductivity,  Btu/ft-sec-°R 
mass  of  the  sleeve  material,  pounds/ft^ 
Nusselt  number,  hx/k 
Prandtl  number,  pCpg/k 

Reynolds  number,  pVx/|i 


recovery  factor, 


T  -  T 

eq  00 

T  -  T 
stag  0 


Stanton  number,  Nu/RPr 


effective  stream  air  temperature  at  the  wall,  some  tempera¬ 
ture  which  gives  a  thermal  potential  which  is  independent 
of  conductance  h,  °R 

stagnation  temperature,  °R 

inside  surface  temperature  of  the  sleeve,  °R 
free-stream  temperature,  °R 
free-stream  velocity,  ft /sec 
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x 
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longitudinal  distance  along  sleeve,  ft 
free-stream  density  of  the  air,  slugs/ft-^ 
dynamic  viscosity  coefficient,  lb-sec/ft^ 

time  derivative  of  average  wall  temperature  for  a  particular 
x-station,  °R/sec 


APPARATUS  AND  METHOD 


The  apparatus  consisted  of  a  three-dimensional  axially  symmetric 
nozzle  which  was  directly  connected  to  the  settling  chamber  of  a  cold-air 
blow-down  jet.  The  nozzle  consists  of  a  shaped  plug  and  an  outer  sleeve 
of  8-inch  extra-heavy  seamless  carbon-steel  pipe  which  was  machined  inside 
and  outside  to  a  wall  thickness  of  O.388  inch.  A  detailed  drawing  of  the 
apparatus  is  shown  in  figure  1.  The  ordinates  of  the  plug  are  given  in 
table  I.  The  plug  was  designed  by  using  three-dimensional  characteristics 
(see  ref .  1) . 

Details  of  the  installation  of  the  thermocouples  and  static-pressure 
orifices  are  shown  in  figures  2  and  3.  The  thermocouples  are  located 
0.060  inch  from  the  inner  surface  of  the  sleeve.  The  wires  are  number- 
30-gage  copper-constantan  (0.010  inch  in  diameter).  For  these  wires, 
conduction  along  the  wire  length  is  negligible.  As  indicated  in  figure  2, 
the  thermocouples  are  in  intimate  contact  with  the  steel  wall  so  that 
thermal  resistance  at  this  junction  is  negligible. 

Figure  1  shows  the  location  of  the  static-pressure  orifices.  The 
Mach  number  distribution  was  measured  both  lengthwise  and  around  the 
sleeve.  The  results  are  shown  in  figure  4. 

The  temperature-recording  equipment  consisted  of  three  synchronized 
high-speed  Brown  recorders,  having  a  total  of  36  switches,  two  of  which 
were  connected  to  the  settling-chamber  thermocouples,  the  rest  being  con¬ 
nected  to  the  sleeve  thermocouples.  The  stagnation-temperature  variation 
with  time  for  a  settling-chamber  pressure  of  220  pounds  per  square  inch 
gage  is  shown  in  figure  5*  Figure  6  shows  the  wall-temperature  variation 
with  time  for  stations  12  and  18  for  a  settling-chamber  pressure  of 
220  pounds  per  square  inch  gage.  Test  runs  were  made  for  settling-chamber 
pressures  of  92,  142,  220,  and  290  pounds  per  square  inch  gage.  Excluding 
the  first  20  seconds,  the  pressures  were  maintained  constant  for  each  run. 
The  recorders  were  calibrated  immediately  before  and  after  each  run  and 
are  accurate  to  *1°  F . 
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In  figure  7  there  are  plotted  for  various  times  during  the  run  the 
curves  of  wall  temperature  against  longitudinal  distance  along  the  cylin¬ 
der.  These  curves  were  used  to  determine  the  rate  of  change  of  the  longi¬ 


tudinal  conduction 


along  the  cylinder. 


only  for  the  length  of  the  cylinder  for  which 


Test  results  were  taken 


was  zero. 


REDUCTION  OF  DATA 


Short  test  runs  in  which  wall-temperature — time  histories  are 
obtained  give  the  essential  data  from  which  heat -transfer  coefficients 
and  recovery  factors  may  be  determined.  The  reduction  of  these  data 
requires  a  method  in  which  the  stagnation  temperature  and,  accordingly, 
the  equilibrium  temperature  may  vary  in  an  arbitrary  manner. 

The  Stanton  number  is  calculated  from 


St 


h 

pVCpg 


and  the  heat-transfer  coefficient  is  calculated  from 


(1) 


h  =  me 


dT/dt 

T  -  T 
xw  xeq 


By  definition,  the  recovery  factor  is 


V  = 


T  -  T 
eq  00 

T~  ^"tI 

stag  “ 


(2) 


(3) 


The  method  of  reducing  the  data  is  simply  to  select  a  recovery  fac¬ 
tor  and  from  equation  (3)  obtain  Te^.  Substitute  this  Te^  into  equa¬ 
tion  (2)  and  obtain  values  of  h  for  different  heat  flow  rates.  These 
values  of  h  are  then  substituted  into  equation  (l).  The  true  values 
of  Tgq,  rjr,  and  St  are  obtained  when  the  St  is  constant  with  time 

(different  heat  flow  rates).  Figure  8  shows  the  curves  used  in  evaluating 
the  Stanton  number  and  the  recovery  factor  for  station  12  for  a  settling- 
chamber  pressure  of  220  pounds  per  square  inch  gage. 
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The  values  of  specific  heat  and  specific  gravity  of  the  sleeve 
material  were  taken  from  reference  2.  Figure  9  gives  the  variation  of 
the  specific  heat  of  the  sleeve  material  with  temperature.  The  values 
of  viscosity  and  Prandtl  number  (0.71)  for  air  were  also  taken  from 
reference  2  and  were  based  upon  the  inside  surface  temperature  of  the 
sleeve.  The  value  of  Tw  chosen  was  at  120  seconds  after  starting. 
(This  arbitrary  choice  of  temperature  can  be  made  because  the  viscosity 
coefficient  is  practically  a  constant  factor  of  the  heat -conduct ion 
coefficient . ) 


ACCURACY 


The  basic  data  consist  of  temperature -time  histories  at  each  thermo¬ 
couple  location.  The  recorders  were  calibrated  for  synchronization  imme¬ 
diately  before  and  after  each  run  and  are  accurate  to  *1°  F.  The  thermo¬ 
couple  wires  are  number-30-gage  ( 0.010- inch  diameter)  copper-constantan 
wires.  Conduction  along  the  length  for  these  wires  is  negligible.  The 
wires  are  insulated  except  at  the  junction  where  intimate  contact  with 
the  sleeve  keeps  thermal  resistance  to  a  minimum. 

The  sleeve  is  enclosed  in  a  vacuum  jacket  to  avoid  free  convection 
and,  for  the  temperature  range  considered,  radiation  effects  range  up  to 
2  percent.  The  product  of  the  mass  per  unit  area  and  the  specific  heat 
f  of  the  sleeve  material  as  used  in  equation  (2)  is  accurate  to  -5  percent. 

The  value  of  dT/dt  used  in  equation  (2)  is  that  of  the  actual 

"  thermocouple  reading,  whereas  it  should  be  the  value  associated  with  the 

average  temperature  in  the  thickness  of  the  sleeve  at  each  point.  During 
the  test  runs,  the  inner-surface  temperatures  are  lower  than  the  average 
temperatures  by  approximately  1°  to  1.5°  at  the  beginning  of  the  runs 
(40  seconds  after  starting)  and  become  the  same  at  the  end  of  the  runs 
(200  seconds).  This  effect  results  in  the  inner-surface  temperature-time 
derivative  being  less  than  the  average  temperature-time  derivative  by 
approximately  1  percent.  The  computed  slopes  when  integrated  reproduce 
the  original  curves  to  within  *0.5°  F. 

The  analysis  of  reference  3  was  used  to  determine  the  effect  of  the 
sleeve  diameter.  This  effect  was  computed  to  be  less  than  0.01  percent. 


RESUITS  AND  DISCUSSION 


Figure  7  shows  the  variation  of  wall  temperature  with  longitudinal 
distance  along  the  cylinder  for  a  settling-chamber  pressure  of  220  pounds 
per  square  inch  gage.  Over  the  test  range,  the  wall  temperatures  are  con¬ 
stant  with  x.  In  the  relation 
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m  and  c  are  constant  and  dT/dt  is  constant  because  Tw  is  constant. 
Therefore,  if  there  is  to  be  a  variation  of  h  with  Reynolds  number  or 
x,  Teq  must  vary  with  x.  The  value  of  Teq  obtained  for  the  run  with 

a  chamber  pressure  of  220  pounds  per  square  inch  gage,  evaluated  at 
120  seconds  after  starting,  actually  decreases  approximately  2°. 


Figure  10  shows  the  variation  of  local  Nusselt  numbers  with  local 
Reynolds  numbers.  The  x  used  in  evaluating  these  numbers  was  based 
upon  an  x  that  was  considered  zero  at  the  minimum.  A  single  line 
faired  through  all  the  data  points  has  a  slope  of  approximately  eight- 
tenths  but,  if  lines  are  faired  through  the  points  for  the  individual 
runs,  the  slopes  are  slightly  higher.  It  is  reasonable  to  assume,  how¬ 
ever,  that  the  higher  stagnation  pressures  would  produce  earlier  transi¬ 
tion  and  hence  the  zero  value  of  x  would  move  upstream  with  increasing 
stagnation  pressure.  Hence,  if  it  is  assumed  that  the  individual  runs 
should  merge  into  a  continuous  line  (rather  than  a  series  of  steps),  x 
may  be  adjusted  for  x  =  0  locations  which  would  make  all  the  runs  coin¬ 
cident.  This  adjustment  has  been  done  in  figure  11  by  using  values  of 

x  equal  to  zero  at  2.00,  4.25,  7-25>  and  10.00  inches  ahead  of  the  mini¬ 
mum  for  the  runs  with  chamber  pressures  of  92,  142,  220,  and  290  pounds 
per  square  inch  gage,  respectively.  This  systematic  adjustment  with 
pressure  of  the  x  =  0  location  does  not  affect  the  value  of  a  faired 
curve  through  all  the  data  (figs.  10  and  11). 

The  Nusselt  numbers  were  found  to  vary  from  4,125  to  32,500  for  the 

Reynolds  number  range  of  5*6  X  I06  to  6.5  X  107  (based  on  adjusted  x  =  0 

locations) .  For  comparison,  the  curves  based  upon  the  analyses  of  Van 
Driest  (ref.  4)  and  Eber  (ref.  5)  are  shown.  In  the  Van  Driest  analysis, 
T-w/Too  was  considered  to  be  3*0  (the  average  value  of  the  test  results). 

For  comparison  with  these  references,  the  data  were  computed  by  using  the 
free-stream  temperature  to  determine  the  density  and  the  velocity.  The 
wall  temperature  was  used  to  determine  the  viscosity  and  Prandtl  number. 
The  two  experimental  values  from  V-2  data  (ref.  6)  are  also  determined 
in  this  manner. 


The  results  of  this  paper  are  approximately  10  percent  higher  than 
the  Van  Driest  analysis  and  approximately  27  percent  higher  than  the 
Eber  analysis.  The  V-2  experimental  results  coincide  with  the  results 
of  this  paper. 

Figure  12  shows  the  variation  of  local  recovery  factors  with  local 
Reynolds  numbers.  The  variation  is  from  O.887  at  5*6  X  10^  to  0.868  at 
6.5  X  107.  The  results  are  compared  with  an  extrapolated  curve  which 
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represents  the  fairing  of  the  data  of  reference  7*  The  data  of  refer¬ 
ence  7  are  for  Reynolds  numbers  below  6.75  X  10^.  The  results  of  this 
report  and  those  of  reference  7  agree  to  within  *0.5  percent.  Also 
included  for  comparison  are  the  curves  for  recovery  factor  equal  to 

Pr^/3  and  Pr"*"/^.  The  wall  temperature  was  used  to  determine  the 
Prandtl  number. 


CONCLUDING  REMARKS 


It  can  be  concluded  that  the  method  of  testing  and  reducing  the 
data  gives  results  that  are  in  good  agreement  with  established  theoretical 
and  experimental  work.  The  reduction  of  the  heat-transfer  data  requires 
adjusting  the  x  =  0  locations  (x  is  the  longitudinal  distance  along 
sleeve)  so  that  they  are  consistent  when  the  Reynolds  numbers  are  varied 
either  by  changing  the  pressure  or  the  "x  value"  location.  The  varia¬ 
tion  of  Nusselts  number  with  Reynolds  number  before  and  after  adjustment 
of  the  x  =  0  locations  are  essentially  the  same. 

The  recovery- fact or  results  are  in  good  agreement  with  the  extra¬ 
polated  results  of  lower  Reynolds  number  data. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  June  24,  1954. 
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Figure  4.-  Mach  number  distribution  for  settling  chamber  pressure  of 

220  pounds  per  square  inch  gage . 
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Figure  5*-  Stagnation  temperature  variation  with  time  for  a  settling 
chamber  pressure  of  220  pounds  per  square  inch  gage. 


Station  12 


Figure  6.-  Wall  temperature  variation  with  time  for  a  settling  chamber 
pressure  of  220  pounds  per  square  inch  gage. 
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Figure  10.-  Variation  of  local  Nusselt  number  with  local  Reynolds  number 

for  x  =  0  at  station  0. 


Nusselt  number 
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Reynolds  number,  R 

Figure  11.-  Variation  of  local  Nusselt  number  with  local  Reynolds  number 
for  corrected  x  =  0  locations.  Viscosity  and  Prandtl  number  determined 
for  wall  temperatures. 


